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Abstract 

For the voter model, we study the effect of a memory-dependent transition rate. 
We assume that the transition of a spin into the opposite state decreases with the 
time it has been in its current state. Counter-intuitively, we find that the time 
to reach a macroscopically ordered state can be accelerated by slowing-down the 
microscopic dynamics in this way. This holds for different network topologies, in- 
cluding fully-connected ones. We find that the ordering dynamics is governed by 
two competing processes which either stabilize the majority or the minority state. 
If the first one dominates, it accelerates the ordering of the system. The conclusions 
of this Letter are not restricted to the voter model, but remain valid to many other 
spin systems as well. 

PACS: 02.50.Ey 64.60.De, 89.65.-s 

How fast an out-of-equilibrium system reaches an ordered state has been a central ques- 
tion in statistical physics, but also also in disciplines such as chemistry, biology, and 
social sciences. Despite its simple structure, the "voter model" has served as a paradigm 
to study this question [l| . It is one of the few spin systems that can be analytically solved 
in regular lattices jl]. In physics, it was investigated how the time to reach the equi- 
librium state depends on the system size, the initial configuration, and the topology of 
the interactions [3|]. Among its prominent properties, the magnetization conservation has 
been studied extensively [4]. Furthermore, the formation and growth of state domains was 
studied, showing the existence of coarsening without surface tension in two-dimensional 
systems {sf. The voter model also found numerous interdisciplinary applications, e.g. in 
chemical kinetics jsl and in ecological and »ciaUy.ten.S. Its properties have also 
served to complete the understanding of other spin-systems, such as the Ising model and 



spin-glasses [10 1 



To assume that transition rates are constant in time is (in general) not valid for non- 
equilibrium systems. A good example are spin glasses, where the effective temperature 



of the system changes with the time elapsed since a given perturbation was applied [H 
In this Letter, we consider that, for each site, the transition rates are not constant, but 
decrease with the time elapsed since the last change of state (namely, its persistence time). 
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We refer to this change as increasing inertia. The level of inertia is measured by how fast 
the transition rates decrease with persistence time. Dependent on the context of the voter 
model, this mechanism has different interpretations. In a social context, the longer a voter 
already stays with its current state, the less it may be inclined to change it in the next 
time step, which can be interpreted as conviction. In models of species competition 0], 
this would imply that neighboring species are less likely to be displaced at a later stage 
of growth. 

Obviously, increasing inertia leads to slower microscopic dynamics. Against intuition and 
in contrast to results with fixed (homogeneous or heterogeneous) values of inertia, we find 
that the time to reach an ordered state can be effectively reduced. We further find that 
this phenomenon exists independently of the exact network topology in which the system 
is embedded. We show that the unexpected reduction of the time to reach an ordered state 
is related to the break of magnetization conservation, which holds for the standard voter 
model. This break originates from the evolving heterogeneity in the transition probabilities 
within the voter population, which, in the extended model, depends on the distribution 
of the persistence times. 

The voter model denotes a simple binary system comprised of voters, each of which 
can be in one of two states (often referred to as opinions), ai = ±1. A voter is selected at 
random and adopts the state of a randomly chosen neighbor. After A^ such update events, 
time is increased by 1. In this work, we consider homogeneous networks, where all voters 
have the same number of neighbors. In the standard voter model, the transition rate at 
which voter i switches to the opposite state, u;^(— (Tj|(Tj), is proportional to the frequency 
of state — (Tj in {i}, the set of the k neighbors of i, namely 



The prefactor (3 determines the time scale of the transitions and is set to /3 = 1. In order 
to describe the dynamics on the macrolevel, we introduce the global densities of voters 
with state +1 as A{t) and with state —1 as B{t). The instantaneous magnetization is 
then given by M{t) = A{t) — B{t). Starting from a random distribution of states, we have 
M(0) = 0. The emergence of a completely ordered state (which is often referred to as 
consensus) , is characterized by \M\ = 1. The time to to reach consensus, T^, is obtained 
through an average over many realizations. The dynamics of the global frequencies is 
formally given by the rate equation 




(1) 



A{t) = -B{t) = - l)B{t) - Q 



{-l\ + l)A{t). 



2m 
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The macroscopic transition rates f2 have to be obtained from the aggregation of the 
microscopic dynamics given by Eq. ([T]). A simple expression for these can be found 
in the mean-field limit. There, it is assumed that the frequencies of states in the lo- 
cal neighborhood can be replaced by the global ones. This gives — 1) = A{t), 
+ 1) = B{t) and leads to A{t) = A{t)B{t) - B{t)A{t) = 0. For an ensemble aver- 
age, the frequency of the outcome of a particular consensus state -|-1 is equal to the initial 
frequency A{0) of state +1, which implies the conservation of magnetization. It is worth 
noticing that, for a single realization, the dynamics of the voter model is a fluctuation 
driven process that, for finite system sizes, always reaches consensus towards either -|-1 
or —1. We now investigate how this dynamics changes if we modify the voter model by 
assuming that voters additionally have an inertia z/j G [0, 1] which leads to a decrease of 
the transition rate to change their state 



Obviously, if all voters have the same fixed value of inertia u,, the dynamics is equivalent 
to the standard voter model with the time scaled by a factor (1 — z/,)~^. Similar results are 
obtained if the inertia values are randomly distributed in the system: higher consensus 
times are found for increasing levels of inertia. In our model, however, we consider an 
individual and evolving inertia z/j that depends on the persistence time the voter has 
been keeping its current state. For the sake of simplicity, the results presented here assume 
that the individual inertia z/j increases linearly with persistence time r^, /i being the 
"strength" of this response, until it reaches a saturation value Ug, i.e. z^(rj) = min [/iTj, Ug]. 
Choosing z^^ < 1 avoids trivial frozen states of the dynamics 0. The rate of inertia growth 
fi determines the number of timesteps until the maximal inertia value is reached, denoted 
as Tg = [us/fA- 

Increasing /i increases the level of inertia within the voter population, thereby slowing- 
down the microscopic dynamics. Like in the case with fixed inertia, one would intuitively 
assume an increase of the average time to reach consensus. Interestingly, this is not always 
the simulation results of Tk(/U.) show for different network topologies (see Fig. [1]). 

Instead, it is found that there is an intermediate value /i*, which leads to a global minimum 
in [^. For yU < /i*, consensus times decrease with increasing values. Only for /i > yU*, 
higher levels of inertia result in increasing consensus times. 



^The results presented here are qualitatively independent of the exact functional relation Vi{Ti), as 
long as a monotonously increasing function with a saturation below 1 is considered. 

^In this Letter, we do not investigate the origin of the global maxima in the consensus times of Fig. [T] 
(a). In contrast to the global minima, this effect results from spatial configurations as can be learned 
from panels (b) and (c) of the same figure. 




(2) 
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Figure 1: (color online). Average consensus times for varying values of the inertia slope 
/i and fixed saturation value Ug = 0.9. Sample sizes vary between 10'^ — 10^ simulation runs. 
Filled, black symbols always indicate the values of a.t fi = 0. (a) 2d regular lattices 
(ki = 4) with system sizes: (o) N = 100, (A) = 400, (□) = 900. The inset shows how 
consensus time scales with system size in regular lattices at /i = fi*: (o) Id, (x) 2d, (<) 
3d, {-k) 4d. (b) Small-world networks obtained by randomly rewiring a 2d regular lattice 
with probability: (o) = 0, (A) Pr = 0.001, (□) Pr = 0.01, (o) p^ = 0.1, (★) Pr = 1. The 
system size is A^ = 900. (c) Fully connected networks (mean field case, ki = N — 1) with 
system sizes: (o) A^ = 100, (□) A^ = 900, (o) A^ = 2500, (★) A^ = 10^. Lines represent 
the numerical solutions of Eqs. (E]), ([6]), (I7j) with the specifications in the text. The inset 
shows the collapse of the simulation curves by scaling /i and as explained in the text. 

For a two-dimensional lattice, shown in Fig. [U^a), we find /i* oc 1/lnA^. Simulations of 
regular lattices in other dimensions show that the non- monotonous effect on the consensus 
times is amplified in higher dimensionality of the system. Being barely noticeable for d = 1, 
the ratio between T^^fi*) and T^i^fi = 0) (i.e. the standard voter model) decreases for d = 3 
and d = A. We further compare the scaling of with system size A^ for the standard 
and the modified voter model. The first one gives for one-dimensional regular lattices 
{d = 1) oc and for two-dimensional regular lattices {d = 2) oc A^logA^. For 
d > 2 the system does not always reach an ordered state in the thermodynamic limit. In 
finite systems, however, one finds ~ A^. In the modified voter model, we instead find 
that Tf^ifi*) scales with system size as a power-law, T^^i^ii*) oc A^" (see inset in Fig. [I](a)); 
where a = 1.99 ± 0.14 for d = 1 (i.e., in agreement with the standard voter model); 
a = 0.98 ± 0.04 for = 2; a = 0.5 ± 0.08 for = 3; and a = 0.3 ± 0.03 for d = 4. For 
fixed values oi fi > fi*, the same scalings apply. 

In order to cope with the network topology, in Fig. [U^b) we plot the dependence of the 
consensus times for small-world networks built with different rewiring probabilities. 
The degree of each node is kept constant by randomly selecting a pair of edges and 
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exchanging their ends with probabihty p 121]. It can be seen that the effect of reduced 



consensus times for intermediate values of /x still exists and is amplified by increasing 
the randomness of the network. This result implies that the spatial extension of the 
system, e.g. in regular lattices, does not play a crucial role in the emergence of this 
phenomenon. This can be confirmed by investigating the case shown in Fig. [11(c), in 
which the neighborhood network is a fully-connected one (the solid lines correspond to 
a theoretical approximation introduced below). The inset shows the results of a scaling 
analysis, exhibiting the collapse of all the curves by applying the scaling relations fi' = 
|/iln(r/iV) - /ill, and = T«/ ln(iV/Oyu', with r] = 1.8(1), /ii = 1.5(1), ^ = 7.5(1). This 
shows that the location of the minimum, as well as T^, scales logarithmically with A^. 

The fact of reaching a final state faster by decelerating the dynamics microscopi cally 



has some resemblances with the "slower-is-faster" effect discovered in panic research |13 
However, the origin of the phenomenon discussed here is quite different, as we can demon- 
strate by the following analytical approach. First, note that voters are fully characterized 
by their current state ±1 and their persistence time r. Thus, we introduce the global 
frequencies aT-(t), b^-it) for subpopulations of voters with state +1,-1 (respectively) and 
persistence time r. Thus, these frequencies satisfy 

A{t) = Y,ar{t), B{t) = J2br{t). (3) 

r r 

Formally, the rate equations for the evolution of these subpopulations in the mean-field 
limit are given by 

r' 

— |^f2(aT-/|aT-) + f2(&^/|aT-) a,-. (4) 

t' 

Due to symmetry, the expressions for br{t) are obtained by consistently exchanging A ^ B 
and a^- ^h^-. 

Note that most of the terms in Eq. (jl]) vanish because for a voter only two transitions are 
possible: (i) it changes its state, thereby resetting its r to zero, or (ii) it keeps its current 
state and increases its persistence time by one. Case (i) is associated with the transition 
rate f2(6o|aT-), that in the mean-field limit reads VL{hQ\ar) = (1 — v{t)) B{t). B{t) is the 
frequency of voters with the opposite state that trigger this transition, while the prefactor 
(1 — z^(t)) is due to the inertia of voters of class a^- to change their state. For case (ii), 
Q{ar+i\ar) = 1 — Q{bo\ar), since no voter can remain in the same subpopulation. I.e., in 
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the mean-field limit, the corresponding transition rates are Q{ar+i\ar) = A(t) + z/(r)i?(t). 
Therefore, if r > 0, Eq. (jll) reduces to 

= A{t) + u{T-l)B{t)\ar-i{t)-ar{t). (5) 
On the other hand, voters with r = evolve as 

«o(^) = ^^lb{,ao\br)br{t) - ao{t) 

T 

= A{t)[B{t) - lB{t)] -ao{t). (6) 

Due to the linear dependence of the transition rates on inertia, the terms involving u can 
be comprised into Isit) and lA{t), namely the average inertia of voters with state —1 and 
+ 1, respectively, i.e. 

Ut) = E ^i-'M) ; lB{t) = u{T)K{t). (7) 

T T 

Expressions (0, El El) and the corresponding ones for subpopulations b^. can be used to 
give an estimate of the time to reach consensus in the mean- field limit. Let us consider an 
initial state ao(t) = A{0) = 1/2 + N-^ and bo{t) = B{0) = l/2-N-\ i.e. voters with state 
-|-1 are in slight majority. By neglecting fluctuations in the frequencies (which drive the 
dynamics in the standard voter model), these equations are iterated until B{t') < N~'^ (i.e. 
for a system size A^, if the frequency of the minority state falls below A^~^, the absorbing 
state is reached). Then, we assume = t'. The full lines in Fig. [H^c) show the results of 
this theoretical approach, exhibiting the minimum and displaying good agreement with 
the simulation results for large values of /i. For low values of /i, fluctuations drive the 
system faster into consensus compared to the deterministic approach. 

Inserting Eqs. (0 ED into the time-derivative of Eq. yields, after some straightforward 
algebra, the time evolution of the global frequencies 

A{t)=lA{t)B{t)-lB{t)A{t). (8) 

Remarkably, the magnetization conservation is now broken because of the influence of the 
evolving inertia in the two possible states. For u{t) = v, (that includes the standard voter 
model, i^, = 0), we regain the magnetization conservation. Interestingly enough, Eq. ([8]) 
implies that the frequency A{t) grows iff. lA{t)/A{t) > iBif) / B{t). 

When the time dependence of the inertia on the persistence time is a linear one, as 
assumed in this Letter, inserting Eqs. (0 El) into Eq. (jTl) we obtain an equation for the 
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time evolution of /a(^) up to first order in fi: 

iA{t) = A{t) lA{t) + fiA\t) - lA{t) + ar). (9) 

Here, ar = J2t>ts contains all subpopulations with maximum inertia. Eqs. (IHl EI) 
correspond to a macroscopic level description of this model. This system of equations has 
a saddle point, A = B = 1/2, Ia = Ib = /^/S + and two stable fixed points, 

one at A = 1, Ia = Vs and another at -B = 1, Ib = Vs- Note that the saddle point 
is close to the initial condition of the simulations. Neglecting fluctuations, the time to 
reach consensus has two main contributions: (i) the time to escape from the saddle point, 
T<,; and (ii) the time to reach the stable fixed point, Tj; namely ~ Tg + Tf. We 
then linearize the system around the fixed points and calculate the largest eigenvalues 
Xg and Xf (for the saddle and the stable fixed points, respectively) as a function of /i. A 
simple argument shows that Tgj ~ In N/\Xsj{fi)\. At the saddle point, we find Xs{n) = 
a/1 + 20/i + 4/1^ — 2yU — l + 0(/i^), which equals to at /i = and monotonously increases 
with fi. For larger values of /i, where the first order term expansion is no longer valid, 
numerical computations show that continues to increase monotonously with fi. This 
means that for larger inertia growth rates /i, the system will escape faster from the saddle 
point, thereby reducing the contribution Ts to the consensus time T^. On the other hand, 
for /i — s> 0, As vanishes and the system leaves the saddle point only due to fluctuations. 

Near the stable fixed points the contribution of to Eq. ([9]) cannot be neglected anymore. 
We then obtain A/^i = —Ug for /i < l — Ug, whilst A/,2 = — 1 for fi > l — Vg. Interestingly, 
both reflect different processes: the eigenvalue A/^i is connected to voters sharing the 
majority state which are, at the level of inertial to adopt the minority one (signalled 
by A/,1 being constant). For fi > l — Us, the largest eigenvalue A/,2 is related to voters with 
the minority state that are, for increasing n, more inertial to adopt the majority state 
(apparent by the decrease in |A/^2|)- 

The contributions Tg and Tf are two competing factors in the dynamics towards consensus. 
Qualitatively, they can be understood as follows: in the beginning of the dynamics, the 
inertia mechanism amplifies any small asymmetry in the initial conditions. While this 
causes faster time to consensus for (small) increasing values of fi, for sufficiently large 
values of inertia growth, another process outweighs the former: the rate of minority voters 
converting to the final consensus state is considerably reduced, too. It is worth mentioning 
that the phenomenon described here is robust against changes in the initial condition: 
starting from Ia = Ib < '^s, if holds for any initial frequencies of opinions. Conversely, 
starting from A = B = 1/2, it holds for any Ia 7^ Ib- 

Summarizing, we investigated the role of microscopically time-dependent transition rates. 
In particular, we consider that the microscopic transition rates decrease with the time 
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elapsed since the last state change of a given site (called inertia). Counterintuitively, we 
find that intermediate inertia values may lead to much lower times to reach the absorbing 
state, i.e. an accelerated dynamics. It is important to emphasize that this final state is 
not an arbitrary one, but most interestingly, it is always the ordered one. The mechanism 
behind this phenomenon is the existence of two competing processes near the initial 
condition and absorbing states. Due to the general analytical approach taken in this 
Letter, we emphasize that this phenomenon is not restricted to the voter model, but is 
expected to appear near the absorbing states of any spin system, whenever the inertia 
mechanism is present. 
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